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Preface 



This volume is based on papers presented at the Workshop on Membrane 
Computing, WMC 2003, which took place in Tarragona, Spain, in the pe- 
riod July 17-July 22, 2003. This was the Fourth Annual Membrane Computing 
Workshop, and the first one held outside Romania. The first three meetings were 
organized in Curtea de Arge§, Romania - they took place in August 2000 (with 
the proceedings published in Lecture Notes in Computer Science , Vol. 2235), 
in August 2001 (with a selection of papers published as a special issue of Fun- 
damenta Informaticae, Vol. 49, Nos. 1-3, 2002), and in August 2002 (with the 
proceedings published in Lecture Notes in Computer Science, Vol. 2597). 

The 2003 workshop was the second workshop of the Molecular Computing 
Network (MolCoNet) funded by the EU Commission in the Fifth Framework 
Program Information Society Technologies (project number IST-200 1-32008). 
The preproceedings of WMC 2003 were published as Technical Report 28/03 of 
the Research Group on Mathematical Linguistics from Rovira i Virgili University, 
Tarragona, and they were available during the workshop. 

The current volume contains only a selection of the papers from the pre- 
proceedings. Moreover, the selected papers have been significantly modified/ 
improved according to the really vivid discussions that took place during the 
workshop — all the selected papers were additionally refereed. The papers in 
the volume cover all the main directions of research in membrane computing, 
ranging from topics in mathematics and theoretical computer science, to (po- 
tential) applications in biology, sorting, ranking, linguistics, and computer gra- 
phics. Several implementations/simulations on computers, computer networks, 
or electronic reconhgurable hardware are also presented. Thus, the volume is 
a faithful illustration of the current state of research in membrane computing 
(a good source of information about membrane computing is the Web page 
http : / /psystems .disco . unimib . it) . 

The workshop was organized by the Research Group on Mathematical Lingui- 
stics from Rovira i Virgili University, Tarragona, Spain, under the auspices of the 
European Molecular Computing Consortium (EMCC) . The program committee 
consisted of Carlos Martin-Vide (Tarragona, Spain), Giancarlo Mauri (Milan, 
Italy), Gheorghe Paun (Bucharest, Romania, and Tarragona, Spain), Grzegorz 
Rozenberg (Leiden, The Netherlands, and Boulder, Colorado, USA), and Arto 
Salomaa (Turku, Finland). 

The editors are indebted to the contributors and to Springer- Verlag for the 
efficient cooperation in the timely production of this volume. 
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The workshop received financial support from a number of sources: MolCoNet 
Project IST-2001-32008 funded by the European Union, Project TIC2002-04220- 
C03-02 funded by the Spanish Ministry of Science and Technology, and the 
Research Group on Mathematical Linguistics of Rovira i Virgili University. 



November 2003 Carlos Martfn-Vide 

Giancarlo Mauri 
Glreorghe Paun 
Grzegorz Rozenberg 
Arto Salomaa 
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Abstract. We propose here a (biologically inspired) model of P sys- 
tem called proton pumping P system that is a special case of evolution- 
communication P system. In cell biology there are transport mechanisms, 
involving protons. We generalize this idea by considering a few differ- 
ent types of protons. A proton pumping P system is, essentially, an 
evolution-communication P system where a special subset of symbol- 
objects (called protons) is used. In such a system we have simple evolu- 
tion rules (classical evolution rules without target indications), symport 
and antiport rules that exchange some objects (among them, possibly, 
other protons) for a proton; taking inspiration from biology, this partic- 
ular type of antiports is often called proton pumping rules. 

We show that, as expected, the new model is universal, using non- 
cooperative rules, symport and antiport rules of weight one, and enough 
types of protons available for the computation. If we decrease the number 
of types of protons to one or two, then the model is at least as powerful 
as ETOL system, provided that (total) weak or strong priority of antiport 
rules over symport and evolution rules are used. 

Finally, we consider some descriptional complexity measures (again, in- 
spired from biology) for the newly introduced model. 



1 Introduction 

In this paper we investigate proton pumping P systems. They are evolution- 
communication P systems, [3], with some restrictions inspired by the biology (in 
what follows we refer to [5] for the elements of membrane computing, to [1] and 
[6] for the elements related to cellular biology). 

* This author’s work was supported by the research grant 2001CAJAL-BURV4 from 
Rovira i Virgili University. 

** This author’s work was supported by the Spanish Ministry of Culture, Education 
and Sport under the Programa National de Formation de Profesorado Universitario 
(FPU) 
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We recall that in the evolution-communication model the computation con- 
sists of two actions: the evolution of the symbol-objects (application of simple 
rewriting rules) and the communication between the regions (application of sym- 
port/antiport rules). 

It has been shown in [2] that evolution-communication P systems are uni- 
versal, using two membranes, non-cooperative evolution rules and symport and 
antiport rules of weight one. 

The proton pumping model is obtained adding a biological restriction to 
the evolution-communication model and in particular over the antiport, rules. 
Considering that, in many bacteria, the only antiports available are those that 
can exchange a proton with some chemical objects, a natural step is to add 
such restrictions to the model. Therefore, a proton pumping P system is an 
evolution-communication P system with a set of special objects called protons 
that are never created and never destroyed and where only antiports that can 
exchange some symbol-objects (also protons among them) for a single proton 
are admitted (inspired by biology, we call such antiport rules proton pumping 
rules) . 

The similarity between the catalyst objects and the protons should be noticed. 
In both cases such special objects are never created and never destroyed; while 
catalysts are used to help some evolution rule to be applied, the protons are used 
to help some communication rules to be applied. 

We show that the proton pumping model is universal, providing that it can 
use sufficient types of different protons during the computation. The proof is 
made, again, using the simulation of programmed grammars with appearance 
checking as in [2]. 

Moreover we show that, when one can use only one or two types of different 
protons (and this is the case in biology) then the proton pumping model is at 
least as powerful as ETOL system, but using (total) weak or strong priority of 
antiport rules over symport and evolution rules. 

Finally, inspired from the fact that in biology the possible different types of 
antiports and symports used are limited (actually there is a constant number of 
them), we have studied some descriptional complexity measures of the systems 
considered. 

2 Definitions 

We start by recalling from [3] the definition of an EC P system. It is given by 
the alphabet, the membrane structure, the multisets of symbol-objects in each 
region, the evolution rules and symport/antiport rules as formalized below. 

Definition 1. An evolution-communication P system (in short, an EC P sys- 
tem,), of degree m > 1, is defined as 



n = (O,p,,w 0 ,wi,w 2 , ■ ■ ■ ■ ■ ■ ,Rm,io), 



where: 
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— O is the alphabet of objects; 

— p is a membrane structure with m membranes injectively labeled 
with 1, 2, • • • , to; 

— Wi are strings which represent multisets over O associated with 

regions 1,2, of p (wq represents the environment); 

— Ri, 1 < i < to, are finite sets of simple evolution rules over O; Ri is asso- 
ciated with the region i of p; a simple evolution rule is of the form u — » v, 
where u and v are strings over the alphabet O; 

— Ri, 1 < i < to, are finite sets of symport/ antiport rules over O; Ri is 
associated with the membrane i of p; 

— i 0 £ {0, 1, 2, • • • , to} is the output region; ifi 0 = 0, then it is the environment, 
otherwise i Q is a label of some membrane of p. 

The basic model assumes that all rules are applied in a nondeterministic, 
maximally parallel way and that there is no priority among the evolution and 
communication rules. The evolutive approach and the communicative approach 
were also proposed, as having (strong) priority of evolution rules and of commu- 
nicative rules, respectively. 

The following notation is used 

NECP m (i,j, a), a £ { ncoo , coo } U {catk | k > 0} 

( PsECP m (i,j , a), a £ {ncoo, coo} U {catk \ k > 0}) 

to denote the family of sets of natural numbers (the family of sets of vectors of 
natural numbers) generated by EC P systems with at most m membranes (as 
usually, to = * if such a number is unbounded), using symport rules of weight 
at most i, antiport rules of weight at most j, and evolution rules that can be 
cooperative (coo), non-cooperative (ncoo), or catalytic (catk), using at most k 
catalysts. 

Now we are ready to give the definition of a proton pumping P system. 
Definition 2. A proton pumping P system of degree to > 1 is defined as 

II (O , P, p, Wq , W\ , • • • , Wm , Rl , * * * , Rm , R\ , * * * , Rjn ' -^1 ) ' ' * > ^o) , (1) 

where (O, p,w 0 ,w 1 , • • • , w m , Ri, ■ ■ ■ , R m , R[ U f?" = Rx, ■ ■ ■ , R! m U R'^ = R m , i Q ) 
is an evolution-communication P system, P C O is the set of protons, R( are 
the sets of symport rules and R " are the sets of antiport, rules (proton pumping 
rules) of the form ( x , out; p, in) or (p, out ; x, in) where x € 0 + and p £ P. 
Every evolution rule is of the form u — > v, where u £ (O — P) + ,v £ (O — P)* . 

The computation of a proton pumping P system evolves like in the case of an 
evolution-communication P system. The m-tuple of multisets of objects present 
at any moment in the regions of II represents the configuration of the system at 
that moment (the TO-tuple (wi, ■ ■ ■ , w m ) is the initial configuration). A transi- 
tion between configurations is governed by the mixed application of the evolution 
rules and of the symport /antiport rules. All objects which can be the “subject” 
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of the rules from the sets R t , 77' , R" , 1 < * < m, 1 < j < m, have to evolve 
by such rules. As usual, the rules from Rj are applied to objects in region i 
and the rules from /?,( and R" govern the communication of objects through 
membrane i. There is no difference between evolution rules and communica- 
tion rules (symports and proton pumping rules): they are chosen and applied in 
the non-deterministic maximally parallel manner. The system continues parallel 
steps until there remain no applicable rules (evolution rules or symport /antiport 
rules) in any region of 77. Then the system halts, and we consider the multi- 
plicities of objects contained in the output region i a , at the moment when the 
system halts, as the result of the computation of 77. The set of all vectors of 
natural numbers computed in this way is denoted by Ps(II). Later we consider 
other ways to define a computation where we introduce some kind of priority 
among the rules involved. 

We use the following notations 

PsProP^(i,j, a), a £ { ncoo , coo} U {catk \ k > 0}, 

to denote the family of sets of vectors of natural numbers) generated by a proton 
pumping P systems with at most m membranes (as usually, m = * if such a 
number is unbounded), k different types of protons (i.e. , k is the cardinality of 
the set P), using symport rules of weight at most i, antiport rules of weight 
at most j, and evolution rules that can be cooperative (coo), non-cooperative 
(ncoo), or catalytic (catk), using at most k catalysts. 

3 Variants: Weak/Strong Priority of (Proton) Pumping 

After we have introduced the basic model of proton pumping P systems, we 
define here two basic variants, derived, to some extent, from biology, and give 
the notions of weak and strong priority as presented in [5] . 

The first variant is with weak priority of proton pumping, where a weak 
priority of proton pumping rules over other kinds of rules is assumed. In this 
case, weak priority means that in the process of assigning rules to objects, first 
the proton pumping rules are assigned in a nondeterministic, maximally parallel 
manner, and then the other rules are assigned to the remaining objects, again 
nondeterministically and in the maximally parallel way. 

This is different from the strong priorities, as usually considered in P systems, 
because the rules with a lower priority, can be also applied in the same step as 
the proton pumping rules, if there are enough objects for them. 

This differs also from the usually considered P systems with priorities, be- 
cause here the priorities are total : they are specified as priorities of (all) proton 
pumping rules over (all) rewriting and symport rules, rather than between indi- 
vidual rules. 

The second variant proposed is with strong priority of pumping. In this case 
the total strong priority is introduced in the following sense: (all) antiport rules 
associated to a membrane have strong priorities over (all) rewriting rules, asso- 
ciated to both regions adjacent to the membrane, and over (all) symport rules, 
moving objects from either region adjacent to the membrane. 
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In other words this means that if a pumping rule is applied in some membrane 
then it blocks all the other evolution and symport rules that could take objects 
from the two regions where the proton pumping rule chooses its objects. 

In case of the weak priority and strong priority variants we use the notation: 

PsProP^ijj, a,wpp), a £ {ncoo, coo} U {catk \ k > 0}, 
[PsProP^i, j, a, spp), a £ {ncoo, coo} U {catk | A: > 0}] 

to denote the family of sets of vectors of natural numbers generated by proton 
pumping P systems with at most m membranes (m = * if such a number is 
unbounded), at most k different types of protons, symport rules of weight at 
most i, antiport rules of weight at most j, and evolution rules that can be 
cooperative (coo), non-cooperative (ncoo), or catalytic (catk) with at most k 
catalysts and weak [strong] pumping priority. 

In what follows we only take in considerations proton pumping P systems 
with non-cooperative evolution rules, in the basic model and in the weak and 
strong priority variants. Moreover, it is known that, in reality, the possible num- 
ber of different types of different antiport, symport and evolution rules used by 
some biological cell is limited by a constant (a fixed number of transport mecha- 
nisms is, for example, available in many bacteria); for this reason we think that it 
is useful to observe some of the biological descriptional complexity parameters of 
the systems considered. In particular, we will be interested in the total number 
of antiport, symport and evolution rules used by a proton pumping system. A 
detailed survey of these types of results obtained in this paper can be found in 
the Appendix. 

Now we need to recall some preliminaries concepts and some notations. First, 
we introduce a useful normal form for ETOL systems. 

Lemma 1. For each L £ ETOL there is an extended tabled Lindenmayer system 
G = (V,T,H,wq) with 2 tables (H = {hi,hi}) generating L, such that the 
terminals are only trivially rewritten: for each a £ T if (a —> a) £ hi U hi, then 
a = a. 

Proof. Let L £ ETOL. Then, there exists an ETOL system Go such that 
L(Gq) = L. For Go, an equivalent ETOL system G\ can be constructed, where 
all terminals are only trivially rewritten. For G\, an equivalent ETOL system G, 
containing only two tables, can be constructed. Moreover, the transformation 
can be performed only on the nonterminals, leaving terminals as they are. □ 

We also give a slightly modified definition of programmed grammars with 
appearance checking, similar to the one used in [4]: 

Definition 3. A programmed grammar (with appearance checking) is a system 
G = (N, T, S, P) , where N is a finite set of nonterminal symbols, T is a finite 
set of terminal symbols (N C\T = ®), S £ N is the start symbol and P is a finite 
set of tuples of the form (r : a —> (3, <j(r), p(r)), where r is a label of a rewriting 
rule a — > (3, 

Lab(P) = {r \ (r : a — > (3, o(r),ip(r))} 
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is the set of labels of rules in P, and a,<p : Lab(P) — > 2 Lab W) ■ cr(r),ip(r) are 
called the success field and the failure field of r, respectively . 

Definition 4. The language generated by a programmed grammar. 

Let (r : a — > (3, a(r),ip(r)) £ P. We say that w' is derived from w in one step by 
applying or skipping the rule with label r (w => r w' ) if either w = xay, w' = xfly 
or w = w' , a Sub(w). In the derivation, pairs of label and word are considered: 
( r,w ) =4> ( r',w ') if w => r w' and either a £ Sub(w) andr' £ tr(r), or a Sub(w) 
and r' £ tp(r). In other words, if a is present in the sentential form, then the 
rule is used and the next rule to be applied is chosen from those with label in 
<r(r), otherwise, the sentential form remains unchanged and we choose the next 
rule from the rules labeled by some element of tp(r). Let =>* be a reflexive and 
transitive closure of=>. The language generated by a programmed grammar G is 
L(G) = {x £ T* | (r, S) =>* (r 1 , w'),w' => r ' x}. 



Remark 1. In this definition it is natural to have w' => r i x rather than (r ; , w') => 
(r",x) because we need not to have the next rule after we have obtained the 
terminal string. If w' = uAv, u,v £ T* , A £ N, (r' : A y,a(r'),(p(r')) £ 
P, and (r,S) =>* (■ r',w then we say that x = uyv belongs to the language 
L(G), even if cr(r') = 0. This definition is family-equivalent to the one with 
(r',w r ) => (r",x) because for any such grammar we could add a dummy rule 
(r : S — > S, 0, 0) to P , and add r to the success and failure fields of all terminal 
rules without changing the language. We take advantage of this fact in the 
universality proof. 

If ip{r) = 0 for each r £ Lcib(P), then the grammar is said to be without 
appearance checking. If a(r) = <p(r) for each r £ Lab(P), then the grammar is 
said to be with unconditional transfer, in such grammars the next rule is chosen 
irrespective of whether the current rule can be effectively used or not). 

Note 1. From now on by programmed grammars we will assume programmed 
grammars with appearance checking with context-free rules. 



Remark 2. In the universality proof the programmed grammars with appearance 
checking will be simulated, considering pairs 

(S = wi,p 0 ), ■ ■ ■ , (w m ,p m -i), ( w m+1 = x,p m ) for 
(.Pi i hi — W 1 ) * * * = r > ( Pmi IXrn) ■> W m UJ m - f_l — X . 

Here, the rule is chosen during the step when it is applied/skipped, rather than 
one step before, po is a new symbol - a starting point of the control sequence. 

During the proof of all the theorems presented in this paper, we will indicate 
with the label ( nx ) the rule (evolution or antiport/symport rule) present in 
equation n and corresponding to the alphabetical place x (for example: (126) 
means “the second rule in equation (12)”). 
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4 Universality of Proton Pumping P Systems 

In this section we give a universality theorem for proton pumping P systems. 
We prove that such systems are universal when they can use sufficient types 
of protons during the computation. The proof is based on the simulation of 
programmed grammars with appearance checking. We recall the following lemma 
(Theorem 1.2.5, [4]): 

Lemma 2. The programmed grammars with appearance checking generate ex- 
actly the family of recursively enumerable languages. 

Finally we recall that the notation Perm(x) indicates the set of all strings that 
can be obtained as a permutation of the string x. 

Theorem 1. PsProPf (1, 1, ncoo) = PsRE. 

Proof. Let L £ RE. Then there exists a programmed grammar with appearance 
checking G = (N, T, P , S) generating L, where Lab{P) = {i | 1 < i < m} and 
N = {Xi | 1 <j i < n}. We will need the following notations: N' = N U {h}, 
N' = (X | X £ X We now define the morphism 7 : (N' U T)* — > (N' U T)* 
by j(x) = x,x£ N r , and 7 ( 2 ;) = x, x € T. Let us construct the P system 

n = (o, p, i-i . s, wi,w 2 , w 3 , Ri,r 2 , R 3 , R’l^R'i- R 3 i 0, R r< 2 > R'l > o), 

P = {Pi, di\i£ Lab(P)} U {r.j \ X ;j £ N} U {p 0 , c, F}, 

N’ = {X | X £ N'}, 

C = {Hf.dP | i £ Lab(P), 1 < j < 5} U {/, g, H}, 

M = {tj I -2 < j < 5} U {#, K, s}, 

A = {ej,e'j | Xj £ N } U {fj | 1 < j < 5}, 

O = PUTUN' \JN' uWuCuAuM, 

h = Ii [213] 3] 2] i> 
wi = Psrir 2 • • • r n , 
w 2 = PocK , 

w 3 = Sh-pip 2 ■ ■■ p m q i <?2 • • • q m , 

Pi = {X ->• 7 (x)df ] | (i : {X x,a(i),ip(i))) £ P} 

U {h hbf^ejfe \ (i : (Xj x,a(i),<p(i))) £ P} 
U{X^#\X£N}U{fj^ fj—i | 2 < j < 5} 

U {bf ] > #, df ] -► # | i£ Lab(P)} U {# ^ #}, (2) 

R 2 ={X ^ Xt 5 \X £ N'} U {tj -£ tj-x | -1 < j < 5} 

U { 6 - 5) b[%,4 6) -+ 4% I i e Lab(P)} 

U {b\ j) -> It; 1 'id ;' 1 -► 4 i_1) | * e Lab(P), 2 < j < 4} 

U (e^ ->■ e' | Xj e X}, 



(3) 
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R 3 = {X -> X \X eN'}\J{g^ H,K -»• K}, (4) 

R[ = {(a, out) | a £ T}, (5) 

R' 2 = {(X, out), (X, in) | X £ N'} U {(e j5 in) \ Xj £ AT}, (6) 

i?' ={(I,m)|IeiV'}U{( 9l in)}, (7) 



R 2 = {( Pj,out ; df\in), ( pj,out; b[ 5 \in) \ i £ o(j)} 

U out; df\in), ( qj,out ; bf ] ,in) \ i £ tp(j)} 

U {(po, out-, df\in), ( p 0 , out', bf\in) \ i £ Lab(P)} 

U {(l, out; pj, in) \ j £ Lab(P) U {0}} 

U {(l, out; qj, in) | j £ Lab(P)} 

U {(e' , out; tj , in) , (rj,out; fi,in) \ X. t £ Nj 

U {(c, out; s, in), ( s , out; F, in)}, (8) 

i?" = {(X, out; c, in) \ X £ N'}, 

U {(pi,out;d^ ,in),{qi,out;bf\in) \ i £ Lab(P)} 

U {(R, out; Pi, in) \ i £ Lab(P) U {0}} 

U {(H, out; qi,in) \ i £ Lab(P)} 

U {(Xj,out;rj,in),(rj,out; K,in) \ Xj £ N} 

U {(c, out; t- 2 , in), (X, out; F, in), ( F , out; I\, in)}. (9) 

The system 77 simulates the programmed grammar G (see also [2]) using the 
following idea. The protons ( pi , qi) serve the role of remembering the label of the 
rule that has been previously applied ( pi ) or skipped ( qi ), appearance checking 
( Tj ), po is the starting point of the control sequence, c is used to sequentialize 
the application of the rules, and F checks that all nonterminals were rewritten 
at the end. 

We also use objects associated with each terminal (T), associated with each 
non-terminal (TV') , including also the rule application failure symbol (h), their 
“first” versions ( N r ), their intermediate versions (TV') , control sequence symbols 
( b\ J 1 if the rule with label j was skipped and d\ R if the rule with label j was 
applied), appearance checking symbols (e^, e'). Objects fj are used to return 
the appearance checking protons, l,g,H are used to return the control sequence 
proton, tj are delaying c for synchronization, K helps F and r,; to block the 
computation, s is used to end the simulation of derivation, switching to checking 
that no more nonterminals are in region 3, and ^ is the trap symbol. 

In region 1 (the outer one) we simulate the application of the context-free 
rules of G , using simple evolution rules (2a), while the rules (8a-8f) enforce the 
“control sequence”, i.e., transitions from an application of (or from the failure to 
apply) a rule of G to the next one. During the simulation, the nonterminals to be 
rewritten are brought into region 1 by rules (9a, 3a, 6a) and the result is returned 
by rules (6b, 7a), except the terminal symbols, ejected into the environment 
by (5). Rules (4a) remove the bars. The failure to apply a rule is simulated 
by (2b). The appearance checking is enforced by the rules (6c, 3g, 8i, 9f). If 
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appearance checking fails (9f was applied), then (9g) is applied and (4c) blocks 
the computation. If the appearance checking succeeds, then ( 8 j) is applied to 
return the proton. Rules ( 8 k, 81) make sure that the computation halts, and 
that when it does, no nonterminals are present in region 1 , because otherwise 
the computation is blocked by (9i, 9j, 4c). 

The following rules are auxiliary. (2c) blocks the computation if (2a) is not 
applied, and this could happen if there is no rule to rewrite nonterminal X. 
Rules (2d) provide a delay, so that if appearance checking is successful, the pro- 
ton can be returned to region 1. Rules (2e, 2f) block the computation if the 
control sequence is invalid, i.e., if the next rewriting rule is not in the success, 
or failure field of the previous rule applied, or skipped, respectively. ( 2 g) helps 
(2c, 2e, 2f) to block. Rules (3b) organize the delay, so that the new nonter- 
minal will be brought for rewriting after the previous one has been processed. 
Rules (3e, 3f) organize the delay so that the proton, corresponding to the rule 
applied/skipped enters region 2 after the proton, corresponding to the previous 
rule applied/skipped, returns to region 3 using rules (3c, 3d, 8 g, 81r, 7b, 4b, 9d, 
9e). Rules (9b, 9c) replace the “current” label of rule by the “previous” label of 
rule. Finally, (9h) returns c to region 2, so that the next rule can be applied or 
skipped, or to finish the computation. 

We will now proceed to a more detailed explanation of why this construction 
works; in what follows we will indicate with R n x the evolution rule present in 
the set R n at (alphabetical) place x and with R n x i ■ ■ ■ Xk one of the rules present 
in the set of evolution rules R n in position xi or X 2 or • • • Xk (the same idea 
is used for the antiport rules, 7?"x, and symport rules, R' n x). Now we pass to 
the description of valid computations. Applying or skipping a rule of G consists 
of 10 transitions of 77. Suppose that at some step the configuration of 77 is 
[i [ 2 [ 3 w Xhirqi- 3 \ 3 pcKt 2 \ 2 FspriTi\ x v, where 

alph(n) U {p, q} = { Pi , qi \ i £ Lab(P)}, 
alph(p) U {rj = {rj \ Xj £ TV}, 

T 3 £({t_ 2 }U{b\ 1 \d^ 1) \i£Lab(P)})*, 

T 2 £{H,f 1 }* , 

Ti £ ({0 U {ej I Xi £ TV})*, 

with the first two unions being disjoint. This corresponds to the simulation of 
derivation of G and having the sentential form in Perm(wXv), where wX £ TV* 
and v £ T*, and the last rule applied or skipped was j if p = pj or p = qj, 
respectively. 

Suppose that we would like to apply the rule i : ( X -* x, a(i), ip(i)) and 
q = Pi . Let x £ Perm(yz), y £ N*,z £ T*. We get the computation from the 
next table: starting from step 0, after 10 transitions, X will be replaced by y, 
v by vz, p by g, T 3 by T^t- 2 d^\ t 2 by t 2 H , ti by t\1. The only variant not 
considered is Up = q , i.e., the last time the rule with label i was applied, and we 
choose to apply it again in case it is in its own success field. Note that p comes 
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to region 3 before q leaves it, so in that case the computation will have no q in 
steps 0-8, no p in step 10, and the claim remains valid. 



Step 


Region 3 


Region 2 


Region 1 


Env 


Rules 


0-10 


whlTT3 


Kt 2 


FspriTi 


V 




0 


Xq 


pc 






R"a 


1 


cq 


px 






R 2 a 


2 


cq 


pXt 5 






R 2 b , R' 2 a 


3 


cq 


pU 


X 




R 2 b , Rid 


4 


cq 


pt 3 


7 {x)df ] 




R 2 b, R' 2 b, R \ , R' 2 ace 


5 


cq 


*27 {y)df ] 


P 


z 


R 2 bd , R3C1 


6 


cqj{y) 


t\df ] lg 


P 


z 


R?bf, l?3a, R 2 gh , R' 3 b 


7 


cqyg 


tod\ 3) P 


l 


z 


R 2 bf, Rib 


8 


cqyH 


t-id[ 2) p 


l 


z 


R 2 bf,R'3de 


9 


cpqy 


t- 2 d[ 1] H 


l 


z 


i?"6c, R" h 


10 


ypt- 2 df ] 


qcH 


l 


z 


like 0 



Again we start with [ 1 [ 2 [ z wXh'KqT3\jpcKT 2 \ 2 FspriTi\ l v. Suppose that we 
would like to skip rule i : (Xi x, cr(i), an d q = Qj ■ We perform: 



Step 


Region 3 


Region 2 


Region 1 


Env 


Rules 


0-10 


wXttt3 


Kt 2 


Fspn 


V 




0 


hq 


pc 


n 




R3 a 


1 


cq 


ph 


Vi 




R 2 a 


2 


cq 


pht 5 


Vi 




R 2 b, R 2 a 


3 


cq 


Pt 4 


hri 




R 2 b , Rib 


4 


cq 


Pt 3 


(tybP ej f 5 n 




R 2 b , R^bc, R 2 bdf, R\d 


5 


cq 


t 2 \h)bf ) e j 


phn 




R 2 bcg, R3CL , Rid 


6 


cq{h) 


tibf’lg 


mh 




i? 3 a, R 2 be , R"gh, R 2 i, Rid 


7 


cqhg 


t 0 bf ] pri 






R36, R 2 be , Rid 


8 


cqhH 


t-ib^pr. 






i?"de, R 2 be, R 2 j 


9 


cpqh 


t-2b^hH 


n 




R'^ch 


10 


hpt- 2 b - 1 ) 


qcHfi 


rilet 




like 0 



So, starting from step 0 , after 10 transitions, p will be replaced by q, T3 by 
T3t-2b[ 1 \ t 2 by t 2 H fi, n by nZe'. The (useless) case p = q is simulated like 
explained in the last paragraph. 

Terminating computation (no nonterminals are present in region 3 , w = s) 
are simulated as follows: 



Step 


Region 3 


Region 2 


Region 1 


Env 


Rules 


0-2 


hqXiTT3 


PKt-2 


prin 


V 




0 




c 


Fs 




R"fc 


1 




sF 


c 




R 2 l 


2 




F 


cs 




NONE 
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Consequently, ^t(£(G9) = Ps{LI). 



□ 



Remark 3. Using exactly the same construction, we can prove that also the 
proton pumping systems with weak priority of pumping are universal. Thus, 
PsProPf(l,l,ncoo,wpp) = PsRE is also true. 

5 Back to Biology: Using One/Two Kinds of Protons 

In this section we prove several results regarding proton pumping P systems 
where a low number of different types of protons is used during the computation. 
Here we show that proton pumping P systems with weak priority pumping and 
strong priority pumping using one or two different types of protons are at least as 
powerful as PsETOL. When proving the results we also calculate the biological 
descriptional complexity parameters discussed before. 

During this section, the following notations will be used. For a given ETOL 
system, G = (V,T, H,wo) in the normal form given in Lemma 1, we set N = 
V — T and N = {a \ a G N}. The homomorphism 7 : TV U T — 1 IV U T is defined 
as 7(a) = a for each a G T and 7 (a) = a for each a G N. 

Similarly, N' = [a' \ a G N}, TV" = {<*" \ a G N}, N^_ = (o'" | a G N}, 
N"" = {a"" | a G N}, as well as N' = {a' \ a G N}, N" = {a" | a G N}, 
N'" = {o'" | a G N}, N"" = {a"" \ aGN}. 



5.1 Weak Priority of Pumping 

We consider first proton pumping P systems with weak priority of pumping. 
Such systems are able to generate at least PsETOL. The results obtained differ 
from each other in the number of different types of proton used (1 or 2), the 
number of membranes and the number of antiport and symport rules used (this 
is one of the biological parameters discussed before). 

Theorem 2. PsETOL C PsProPl(l,l,ncoo,wpp). 

Proof. Let L G ETOL. Then there exists an ETOL system G = 
(V. T, {hi, /12}, u>o) in the normal form from Lemma 1 generating L. Let us con- 
struct the P system 



n = ( O , {p}, p, s,w ly W2,w 3 ,W4, Rx, R 2 , R 3 , R 4 , R[, R' 2 , R' 3 , R' 4 ,[ 



,R’f,R’f,R'i,0), 



0 = TUNUN'\J N" U TV U {S,p}, 

h = [ 1 [ 2 [ 3 [ 4] 4] 3] 2] 1’ 

Wi = S, 

W 2 = £, 



W3 = P, 



where: 




12 



A. Alhazov and M. Cavaliere 



W4 = e, 

R\ — {A — y ex, ci — )■ y(o) | (A — y cx) £ h\\ 

U {a — y a', a 1 —y a" \ a £ N} U {S — y u>o, S — y 7(^0)}? (10) 

R 2 = R 3 = {a" -> a" | a £ N}, (11) 

R 4 — {A — y cx, a — ^ 7(a) | (A — y ex) £ h 2 } 

U {a — y a' , a' —> a” \ a £ TV}, (12) 

R\ = {(a, out ) | a € T}, (13) 

R' 2 = {(a, out) | a £ T} U {(a, in), (a, out) \ a € N}, (14) 

i ?3 = {(a, out) | a £ T} U {(a, in), (a, out) | a £ N}, (15) 

i ?4 = {(a, out) | a £ T} U {(a, in), (a, out) \ a € N}, (16) 

R '2 ~ {(Pj ou a " ) * n ) I a £ (17) 

i ?3 = {{p, out', a, in), (a, out\ p, in) | a € N}, (18) 

i ?4 = {{a" , out', p, in) \ a £ A^}. (19) 



The idea of the proof is the following. The rewriting rules in region 1 simulate 
the table hi by (10a, 10b), and the rules in region 4 simulate the table h 2 by 
(12a, 12b). The application of rules (10b) : a " —> 7(a) in region 1 or of rules 
(12a) : a" — > a in region 4 will lead to changing the table, as the produced 
symbols will be moved to region 4 by (14b, 15b, 16b), or to region 1 by (16c, 
15c, 14c), respectively. Also, the protons move to the adjacent region because of 
the weak priority of the proton pumping rules. 

To simulate the derivation, all symbols should change the table at the same 
time. We will now explain how the system verifies it after applying the first table 
(the sentential form is in region 1, the proton is in region 3). Let us consider the 
case when both rules (10a) : a" — > 7(a) and (10b) : a” — > a are applied at some 
step s. Then during step s + 1 at least one of the objects, say b, will go into 
the region 2 by (14b) and at least one of the objects will be rewritten to a' by 
(10c). During step s + 2 the proton will be exchanged with b by (18b), because 
of the priority of proton pumping over symport (15c), and a' will be rewritten 
into a" by (lOd). During step s + 3 the proton will be exchanged with a" by (17) 
because antiport has priority over evolution rule, and the computation will be 
blocked by the rule (11) : a" — > a" . 

The case of simulating the application of the second table is symmetric (the 
sentential form is in region 4, the proton is in region 2). The terminal symbols 
are ejected into the environment by rules (13, 14a, 15a, 16a). □ 



Note 2. In the system 77 considered above, the number of different antiport rules 
needed is 41771, where I TV I is the number of nonterminal symbols of the starting 
ET0L system. 



Theorem 3. PsETOL C PsProP^(l,l,ncoo,wpp). 
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Proof. Let L £ ETOL. Then there exists an ETOL system G = 
(V, T, {h\, /i 2 }, Wq) in the normal form from lemma 1 generating L. Let us con- 
struct the P system 

n = (O, {p,q}, {!,,£, W!,W2,W 3 , R,!, R 2 , R3, R[, R'l, R'2, 

where: 

O = T U TV U N' U N" U N U {S,p, q}, 

h = [ 1 [2 [ 3 ] 3 ] 2 ] 1 ’ 
wi = Sq, 
w 2 = £, 
w 3 =p, 

R\ — {n — )■ ex, o — ^ 7 (o) | (^4 — y ex) £ hi} 



U {a — > a , a — > a” \ a £ N} U {S — > w 0 , S — > 7 (^ 0 )}, (20) 

R ‘2 = {a" — >■ a" | a £ N}, (21) 

i ?3 = {n — ^ o, ci — ^ 7 (a) j (A — ^ o) £ h 2 } 

U {a — > a' , a' — > a” \ a £ Nj, (22) 

R'i = {( a J ou t) I a ^ T 1 }) (23) 

R' 2 = {(a, out) | a £ T} U {(a, in), (a, out) \ a £ N} U {(< 7 , out)}, (24) 

i?3 = {(a, out) | a £ T} U {(a, in), (a, out) \ a £ N} U {( p , in)}, (25) 

i ?2 = {(p, out ; a" , in), (a, out ; q, in) | a £ A^}, (26) 

i?3 = {(p, out ; a, in), ( a ", out ; q, in) \ a £ N}. (27) 



The rewriting rules in region 1 simulate hi , and the rules in region 3 simulate 
ft- 2 - The application of rules a" —> 7 (a) in region 1 or of rules a" —> a in region 3 
will lead to changing the table, as the produced symbols will be moved to region 
3, or to region 1, respectively (using the symport rules (24b) and (25b) or (24c) 
and (25c)). 

To simulate the derivation, all symbols should change the table at the same 
time. Let us consider the case when both rules a" — > 7 (a) and a" —> a are 
applied at some step s. Then during step s + 1 at least one of the objects, say 
b, will go into the middle membrane and at least one of the objects will be 
rewritten to a' . During step s + 2 the proton will be exchanged with b and a' 
will be rewritten into a" . During step s + 3 the proton will be exchanged with a" 
because antiport rules have priority over evolution rules, and the computation 
will be blocked by the rule a" — > a" . □ 



Note 3. In the system II considered above, the number of different evolution 
rules needed is 3|1V| + 2r + 3, where |1V| is the number of nonterminal symbols 
and r = |fti| + I/ 12 I is the number of the rules in the ETOL system we started 
with. 
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Theorem 4. PsETOL C PsProPf(l, 1, ncoo, wpp). 

Proof. Let L € ETOL. We start again from an ETOL system G = (V, T, {hi, h 2 }, 
wq) as in Lemma 1 and we construct the P system 

n = (O, {p, q}, p, s, Wi,w 2 , R\,R 2 , 0 , R'l- 0 , R 2 , 1 ), 



where: 



0 = TL)NUN'U N” U N"' U N U TV' U TV" U TV"' U {S, p, q}, 
h — [ 1 [2] 2] i’ 

w\ = S , 
w 2 = pqsk , 

Ri = { at " — > a , a'" — > 7 (a) | {A — > a) € hi} 

U {a'" — > a , a'" — > 7(a) | {A — > a) G h 2 } 

U {a — > a' , a' — > a " , a" a'" | a £ N U N} 



U {/S' — >■ wo, S j(wo), k — > k}, (28) 

R 2 = { a -> a'",a^a 777 \ a <E N}, (29) 

i ?2 = { {a'", out), (a"', out) \ a £ T} U {(s,*?r)}, (30) 

f ?2 = {( p , out ; a, in), ( q , out ; a, in) \ a € N} 

U {(s, out-p, in), (s, out ; q, in), ( k , out;p, in), ( k , out ; q, in)}. (31) 



The rewriting rules in region 1 simulate hi and h 2 . The application of rules 
a 1 " —> 7(a) or of rules a"' —> a will lead to changing the table. 

To simulate the derivation, all symbols should change the table at the same 
time. Let us consider the case when both rules a'" —> 7(a) and a'” —> a are 
applied at some step s. Then at step s + 1 both p and q will be brought in the 
outer region (using the antiport rules (30a) and (30b)) and at step s + 2 both 
s and k will have to come into the outer region (using the rules (31c), (31d), 
(31e), (31f )) in exchange for p and q, and the computation will be blocked by 
the rule k — > k in region 1. 

If only one of the rules is applied, then only one of p and q will be brought in 
the outer region. The symbol a or a will be rewritten to a'" or a , respectively, 
and returns to region 1. The corresponding proton p or q will be exchanged for 
s (if it is exchanged for k, the computation is blocked) and then s returns to 
region 2, and the simulation continues. The computation halts when no more 
nonterminals remain in the system, leaving the result in region 1. □ 



Note f. In the system 17 considered above, the number of different antiport 
rules needed is 2|7V| +4, where \N\ is the number of nonterminal symbols of the 
starting ETOL system. 



Theorem 5. PsETOL C PsProPf (0,1, ncoo, wpp). 
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Proof. Let L G ETOL and an ETOL system G = (V, T, {hi, h 2 }, wo) as in Lemma 
1 generating L. Let us construct the P system 







n 


= (0,{p},p,£,W 1 , 


,W 2 ,Rl,R2, 0 , 0 , 0 ,i? 2 ,l), 




where: 












O 


= 


TUN UN' 


U N" U N"' 


U IV"" 


U TV U TV' U IV" U AT"' U IV"" U 


{S,p,q}, 


h 


= 


[ l [ 2 ] 2 ] i’ 










Wi 


= 


s , 










U>2 


= 


P, 










Ri 


= 


{a"" ->• a, a 


"" -»• 7(a) 1 




0 ) € /ii} 






u 


{a"" — >■ a, a 


"" j(a) | 




a) g Ii 2 } 






u 


{a — > a! , a' - 


a", a" ->• 


a"', a'" 


a"" \ a£ NUN} 






u 


{S w 0 ,S 


-> 7(^o)}, 






(32) 


R 2 


= 


{a a , a' 


— » a 77 ", a 7 " - 


-Go 777 | 


a G IV}, 


(33) 


R'l 


= 


{(p, out ; a, in), ( p , out; a' 


, m) , (p, out; a 1 " , in) , 


(34) 






(a'", out; p, in), (a"", out 


p, in) 


| a G IV}. 


(35) 



The rewriting rules in region 1 simulate hi and /i 2 . The application of rules 
a"" —> 7 (a) or of rules a"" —> a will lead to changing the table. 

To simulate the derivation, all symbols should change the table at the same 
time. 

If only nonbarred letters are present in region 1, then 

[ l wa[ 2 p] 2 ] 1 =>(32e,35 a) [i w, p[ 2 a ] 2 1 1 == ^(32/, 356) \i w "p[2 a '"\ 2 ] i =^(32g,35d) 

[ l w"' a'"[ 2 p] 2 ] 1 =>( 32 / () 35 c) [ 1 'u/" , a""[ 2 p] 2 ] ! and then the first table is applied by 
(32a, 35b). 

If only barred letters are present in region 1, then 

[l^[ 2 p] ah =>(32e) [i^[ 2 p] 2 ]i =>(3 2/, 35b) [ 2 ] i =>(32g,33b) 

iiV'"p{ 2 b""] 2 ] i =>( 3 2 ?i, 35 d) [iU ,,,, 6 ""[ 2 p] 2 ] i and then the second table is applied 
by (32c, 35d). 

If both barred and nonbarred letters are present in region 1, then 
[ 1 wv ab[ 2 p \ 2 \ 1 =>( 32 e, 35 a) [ ± w'v'b'p^ a] 2 ] x =>( 32 /, 33 a) 

[ 1 W , V , 6 >[ 2 a " , ] 2 ] 1 =>(32g,35d) [ l w'" v'"b’" a"' [ 2 p\ 2 ] 1 =>(32M5c) 
[ 1 u;""u""a""p[ 2 ^ ,,/ ] 2 ] ! and then the computation is blocked by the rule (33c). 

The computation halts when no more nonterminals remain in the system, 
leaving the result in region 1 . □ 

5.2 Strong Priority of Pumping 

In this section we show that also proton pumping systems with strong priority 
of pumping are able to generate at least PsETOL , using one type of protons and 
two membranes. 
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Theorem 6. PsETOL C PsProPf(0, 1, ncoo, spp) . 

Proof. Consider an ETOL system G = (V,T,{hi,h 2 },Wo) in the normal form 
from Lemma 1. We construct the P system 

n = (0,{p},fJ,,£,Wi,W 2 ,R 1 ,R 2 ,<I),<I),<!),R^,l), 

where: 

0 = Tl)NUN'\J N" U N'" U N"" 

UNUN 7 LiW 7 uW 77 uW 777 U{S,p,q}, 

h — [ i [ 2 ] 2 ] i> 

wi = S, 

w 2 =p, 

Ri = {a'"' -A a, a'"' — » 7(a) | (A — > a) G hi} 

U {a”" -A a, a"" — > 7(a) | {A — > a) G h 2 } 

U {a -a a', a' -a a", a" -a a'", a'" -a a"" | a £ N U N} 





U {S Wo, S 7 (^ 0 )}, 




(36) 


R‘2 


= { a — y cl , cl — y cl , 








a' -A a", a" -A a, a -A a"", a'" -A a'" a G N}, 


(37) 


to ^ 


= {( p , out ; a, in), (p, out; a' 


, in) , (p, out ; a'" , in) , 


(38) 




(a w , out ; p, in), (a"", out; 


p, in) | a G iV}. 


(39) 



The rewriting rules in region 1 simulate h\ and h 2 . The application of rules 
a "" ^(a) or of rules a"" -A a will lead to changing the table. To simulate the 

derivation, all symbols should change the table at the same time. 

If only nonbarred letters are present in region 1, then 
[ 1 wa[ 2 p\ 2 ] 1 =>(39 a) [ 1 wp[ 2 a\ 2 \ 1 =>( 3 6e,37 a) [ iW'p[ 2 a '] 2 ] 1 =>(36/, 37ft) 
[ 1 w"p[ 2 a"} 2 ] 1 =>(36 d) [i w"a"[ 2 p\ 2 \ 1 =>(36g) [i w"’ a"'[ 2 p] 2 ] 1 =>(36 h) 
l^w"" a""[ 2 p] 2 ] 1 and then the first table is applied. 

If only barred letters are present in region 1, then 

[ApIJi =^(3 9a) [ 1 'V'b%p] 2 ]_ 1 =>(36e,37 a) 2 ] i =>(36/, 376) 

[ 1 u " p |^/ , ] 2 ] 1 =>( 3 9 d ) [ 1 V , " P [ 2 6] 2 ] 1 => (36 e) [ 1 U /, " p [ 2 6""] 2 ] 1 => (39c) 

[ 1 u""6""[ 2 p] 2 ] x and then the second table is applied. 



If both barred and nonbarred letters are present in region 1, then 
[ 1 wvab[ 2 p] 2 ] 1 => (39o) [ 1 wvbp[ 2 a] 2 ] 1 => (36ei 37a) hw'VFp^a'} 2 ] i => ( 36/.37&) 

"v"b"a"[ 2 p\ 2 ] i 






J'l 1 
2 J 2l 1 



> (39d) 



>(36e) 



[ 1 u/ , V"6 , "a" , [ 2 p] 2 ] 1 =>( 3 Q C ) [ 1 u/ , V"a ,,, p[ 2 6 ,,, ] 2 ] 1 and then the computation is 
blocked by the rule (37f): b'" — > V". 

The computation halts when no more nonterminals remain in the system, 
leaving the result in region 1. □ 
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6 Conclusions and Open Problems 

Inspired from the functioning of transport mechanisms and from the proton 
pumping process well-known in cellular biology we have introduced and studied 
a variant of the evolution-communication P systems presented in [3], called 
proton pumping P systems. In this variant we impose a biological restriction 
over the type of antiport rules used by the system: the antiports (here called 
also proton pumping rules) can only exchange some chemical objects (symbol 
objects) for a special object called proton. In a proton pumping P system we 
can have different types of protons (each one with its own number of copies) and 
such special objects are never created and never destroyed; in other words, they 
are only used to activate the antiport rules. 

We have shown how, with enough types of different protons, the proton 
pumping model is universal. Moreover we have introduced two variants of the 
basic model: one with weak priority proton pumping (where the proton pumping 
rules have weak priority over the other rules) and one with strong priority proton 
pumping (where the proton pumping rules have strong priority over the other 
rules). We have shown that these two variants can generate at least the Pariklr 
sets of ETOL languages, using 1 or 2 types of protons. Finally we have considered 
some clescriptional complexity measures such as the different number of types of 
antiport and symport rules used by a system. 

Many interesting problems have been left open. For instance, we know that 
the basic model of proton pumping P system is universal, but with a number 
of different types of proton that is unbounded. Can this number be bounded? 
Moreover, we do not have an upper bound for the two variants considered: they 
can generate at least the PsETOL but we do not know, for example, if they are 
universal. What about considering a proton pumping system with few types of 
protons and few types of catalysts (both are biologically motivated)? 

Finally we notice that in every proof the number of different antiport and 
symport rules used is not bounded and this can be considered not very realistic; 
what about a system with a fixed (a constant) number of different antiports 
or symports available? (What can we generate with such (biological) restricted 
proton pumping systems?) 



Acknowledgements. The authors thank I. Ardelean for the inspiring discus- 
sions. 



References 

1. B. Alberts, Essential Cell Biology, An Introduction to the Molecular Biology of the 
Cell, Garland Publ, New York, London, 1998. 

2. A. Alhazov, Minimizing Evolution-Communication P Systems and EC P Au- 
tomata, New Generation Computing, accepted. 




18 



A. Alhazov and M. Cavaliere 



3. M. Cavaliere, Evolution-Communication P Systems, Membrane Computing (Gh. 
Paun, G. Rozenberg, A. Salomaa, C. Zandron, eds.), LNCS 2597, Springer- Verlag, 
Berlin, 2003, 134-145. 

4. J. Dassow, Gh. Paun, Regulated Rewriting in Formal Language Theory , Springer- 
Verlag, Berlin, Heidelberg, 1989. 

5. Gh. Paun, Membrane Computing. An Introduction, Springer- Verlag, Berlin, Hei- 
delberg, 2002. 

6. M.H. Saier, jr., A Functional-Phylogenetic Classification System for Transmem- 
brane Solute Transporters, Microbiology and Molecular Biology Reviews, 2000, 
354-411. 

7 Appendix 

Here is the summary of the results proved in theorems. Also the slowdown of the 
P systems considered is summarized (i.e., how many steps the P system used 
needs to realize a single step of the simulated grammar) . 



Theorem 


Membranes 


Protons 


Priority 


Class 


Slowdown 


5.1 


4 


1 


wpp 


ET0L 


3-6 


5.2 


3 


2 


wpp 


ET0L 


3 


5.3 


2 


2 


wpp 


ET0L 


4 


5.4 


2 


1 


wpp 


ET0L 


5 


5.5 


2 


1 


spp 


ET0L 


7 


4.1 


3 


2r + |1V| + 3 




RE 


10 




3 


2r + \N\ + 3 


wpp 


RE 


10 



Here is the table of descriptional complexity (biologically motivated) of the 
systems considered: N is the number of nonterminals, T is the number of ter- 
minals, and r indicates the number of productions used by the grammar that 
has been simulated (i.e., r = card(h\) + card(h, 2 ) if the grammar is ET0L, and 
r = card(Lab(P)) in case of programmed grammar). 



Thm 


Antiports 


Uniports 


Evolution 


Objects 


5.1 


4|iV| 


21V + \T\ 


3 1V + 2r + 3 


4|1V| + \T\ +2 


5.2 


4|JV| 


2\N\ + \T\ +2 


3 \N +2r + 3 


4|1V| + \T\ +3 


5.3 


2|iV| +4 


2|1V| + 1 


6|-ZV| + 2r + 3 


8 IV + T + 4 


5.4 


5\N\ 


0 


11 IV + 2r + 3 


10|lV| + |T| + 2 


5.1 


5|iV| 


0 


14 1V +2r + 3 


10 IV j+ T +2 


4.1 


8r + r 2 + 5 IV + 6 


4|1V| + \T\ +4 


|1V| +4r + 8 


12r + 7|lV| + |T|+21 
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Abstract. This paper defines membrane processors as digital processors 
capable of implementing local processing performed inside membranes 
in transition P systems. In order to encode the membrane structures of 
transition P systems, additional binary data structures are needed. Such 
binary data structures are named ’’connectivity arrays”. The main pur- 
poses of the connectivity arrays are to distribute information about the 
system membrane structure among the different processors that imple- 
ment the transition P system, and to facilitate communication processes 
among different membrane processors. The information kept in processor 
has to be versatile and compact; versatile means to easily permit changes 
in the membrane structure of the system processors, and compact means 
that not too much memory space is needed. 

Considering these two requests, a binary word of variable length is con- 
sidered here in order to represent in membrane processors the local in- 
formation about the membrane structures of transition P systems. This 
paper presents the abstract data type “connectivity array” which defines 
the data structure and the needed operations over it in order to parse 
the P system into a membrane processor tree and to keep congruent the 
P system membrane structure with the membrane processor tree. 

1 Introduction 

P systems were introduced by Glreorghe Paun in 1998 [6], and since then this 
research field has experimented an impressive growth. Many theoretical papers 
have been published in different workshops, congresses and scientific journals. 

* Work partially supported by contribution of EU commission under the Fifth Frame- 
work Programme, project “MolCoNet” IST-2001-32008. 



C. Martm-Vide et al. (Eds.): WMC 2003, LNCS 2933, pp. 19—30, 2004. 
(c) Springer- Verlag Berlin Heidelberg 2004 




20 



F. Arroyo et al. 



P systems are computation models which are inspired by basic features of 
biological membranes. Membranes define compartments or regions. Multisets of 
objects and evolution rules are then placed inside regions. Objects evolve by 
means of evolution rules which are applied in a maximally parallel nondeter- 
ministic manner. Objects can pass through membranes, and membranes can be 
dissolved. These are the main features used in P systems. Many variants of P 
systems have been investigated and most of them are computationally universal. 

The information about the membrane structure is very important for imple- 
menting transition P systems, for instance, in order to select evolution rules that 
can be applied in each evolution step and for circulating objects among mem- 
brane processors. Therefore, it is important to keep congruent the membrane 
structure information in membrane processors. 

The next section is devoted to specify the way in which P systems compute. 



1.1 Transition P Systems as Computational Devices 

P systems compute starting from an initial configuration and passing from one 
configuration to another by applying evolution rules to objects inside regions 
until reaching a halting configuration, where there is no rule applicable to the 
objects present in the system. 

A configuration is defined as a (fc+l)-tuple (/z,^, where /z is the 

membrane structure and uj ij is the multiset of objects associated to the region 
defined by the membrane with label ij in the system [7]. 

Given two configurations Ci = (//, ..., uj ' ik ) and C 2 = (/x",w", ..., w") of 

a P system 77 it can be said that there is a transition from C\ to C 2 and this is 
represented by Ci =>■ C 2 , if we can pass from C\ to C 2 by using the evolution 
rules from R tl , ..., Ri k in the regions ii, ..., ik- 

From this definition of computation in P systems, two different levels of 
processing can be defined: 

1. The global process performed by the P system. This process is got by tran- 
sitions between two consecutive configurations. 

2. The local process performed at the membrane level. This process is respon- 
sible for making evolve each region of the system by application of evolution 
rules. 

The global processing is achieved from the local processing by communication 
and synchronization among the different membranes of the system. 

Local processing is performed inside regions, which are defined by the mem- 
branes of the system. Specifically, for every region in the system, the following 
steps are performed in order to achieve the local processing [1,5]: 

1. The computation of a complete multiset of evolution rules. This will permit 
to evolve the region by application of the multiset of rules and will determine 
the multiset of objects that will be sent to the adjacent regions and the 
resulting multiset of the region. 

2. Synchronization and objects communication among adjacent regions. 
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3. Rearranging the membrane structure, as some of the membranes can be 
dissolved. This means that synchronization and communication are needed 
again in order to compute the new system configuration. 



2 Membrane Processors for Transition P Systems 

A membrane processor can be defined as a digital processor able to implement the 
membrane local processing performed inside transition P systems membranes. 

In order to implement transition P systems execution in membrane proces- 
sors, the initial transition P systems membrane structure must be projected to 
the appropriate interconnection of several membrane processors. This projection 
defines a membrane processor tree. 

Moreover, due to the dissolving capability of some evolution rules, the mem- 
brane structure of the transition P system may be changed during the execution. 
However, membrane processors cannot be easily removed from the hardware ar- 
chitecture. One possible solution to this issue is to define two different operating 
modes for membrane processors: 

— The active mode, when the membrane is alive in the system. 

— The passive mode, when the membrane is not alive in the system. 

When a membrane processor is operating in its active mode, it makes evolve 
the system to a new configuration. The following processes have to be performed 
by the processor: 

1. The computation of a complete multiset of evolution rules for each evolution 
step, taking into account the rules and objects stored in the processor. The 
computation is performed according to the following equations: 

labels = {i g L\((i, _, _),_)€ nn} (1) 

ACTIVES = Active(Applicable((Useful R)LABELS)io)p 
COMPLETES = ( Multicomplete ACTIVES)lu 

For a complete description of these equations we refer to [3,5]. 

2. Synchronization with others membrane processors. The processor cannot 
send objects to its adjacent processors if they have not concluded step 1. 

3. Objects communication. Once all processors in the system have finished the 
computation of the complete multiset of evolution rules, the processor has to 
send objects to its adjacent processors. This process produces objects which 
pass through membranes and new objects come to the processor. 

4. Computation of the new multiset of objects. The communication and appli- 
cation of the complete multiset of evolution rules produce a new multiset of 
objects in the processor. 

5. If a membrane dissolution takes place, then the processor has to pass to 
the passive mode; before passing to this mode some processes have to be 
performed: 
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a) Synchronization with the children processors. Due to the fact that some 
of the children processors can dissolve their membranes too, before send- 
ing any information to the father processor, the membrane processor 
must wait for information coming from its children processors. 

b) Computation of the resulting multiset of objects. Once again, a new 
multiset of objects has to be calculated if some of the children processors 
have changed to the passive mode. 

c) Communication with the father membrane processor. At the end of the 
active mode, the processor has to send to the father processor the needed 
information in order to rearrange the system membrane structure. This 
information is related to objects and connectivity as follows:. 

i. Objects multiset. If the membrane dissolution is accomplished, then 
the rules of the membrane are removed from the system and the 
multiset of objects is added to the multiset of object of the father 
membrane processor. 

ii. Connectivity array. The information related to the membrane struc- 
ture is stored into binary words named “connectivity arrays” . Mem- 
brane processors use connectivity arrays for storing different data 
used to rearrange the membrane structure during computation. 

When a membrane processor is computing in the passive mode, it serves as 
communication bridge among others membrane processors. The main function- 
ality of the processor operating in this mode is to send objects and connectivity 
arrays from its children processors to its father processor and viceversa. 

Until now, it has been described what the membrane processors are and their 
different operating modes. The next section is devoted to defining the abstract 
data type named “connectivity array” and to describe how it can be used to 
imply changes of the P system membrane structure to the associated membrane 
processors tree. 



3 Connectivity Arrays 

Connectivity arrays have to accomplish two important tasks. The first one is 
performed while membrane processors are operating in the active mode, and the 
second one is performed when membrane processors are operating in the passive 
mode. 

Every membrane processor, involved in the hardware implementation of a 
transition P system 77 of degree n, stores three different connectivity arrays or 
binary words of length n+1: Processor Identification ( Id-Proc ), Processor Mask 
( Mask-Proc ) and the membrane structure information for the processor ( I 

Let us define more exactly these three binary words. 

Let 77 be a transition P system of Degree(II) = n and let L = {0, 1, ...,n} 
be the set of membrane labels of 77 including the label 0 for the output of the 
P system; then for all i £ L — {0} the following binary words are defined: 
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1. Processor Identification: for all j £ L, 

Id Proc(i)(j) = ^ ^eTwk 

2. Processor Mask: for all j £ L, 

MaskJProc(i)(j) = 4 ^ J < 

[0 otherwise. 

3. 'f'(i) connectivity array: for all j £ L , 

{ 1 if (j is a child of t) V (i = j) V (j is the father of «)V 
(j = 1 A j = 0), 

0 otherwise. 

For example, a typical connectivity array for a membrane processor i indi- 
cating the position of the bits can be represented as follows: 

0 1 2 ri 

= (bi o, b h , b i2 , ■ ■ ■ , 1, ■ • • , b in ) 

Neither Identification Processor nor Mask Processor change during the com- 
putation. They are auxiliary binary words necessary in order to keep the tf'(i) 
connectivity array updated with the appropriate information about the mem- 
brane structure related to the membrane processor i. 

The \P( i ) connectivity array provides to the membrane processor all the in- 
formation about the membrane structure that it needs. However, one restriction 
must be imposed to the order of numbering membranes in the transition P sys- 
tem. The numbering starts with label 1 for the skin membrane and proceeds 
from top down in the membrane processor tree (see Fig. 1) and it provides to 
the active processors information about: 

— their father processor, which is identified by the position of the leftmost bit 
set to one in the connectivity array; 

— the processor itself identified by the second left bit set to one in the connec- 
tivity array while the membrane processor is in the active operating mode; 
— their children processors, which are identified by the position of the rest of 
the bits set to one in the connectivity array. 

During system operation, the if - membrane processor connectivity array may 
be dynamically changed, according to the changes in membrane processors op- 
erating mode. When a membrane processor changes its operating mode from 
active to passive, then it has to change its S' connectivity array and to commu- 
nicate this fact to its father membrane processor and to its children membrane 
processors in order to rearrange their S' connectivity arrays. The way to do this 
is shown below. 

We define several operators in order to properly change the connectivity 
arrays in membrane processors; they are given below in the precedence order: 
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1. <<<: This operator changes to zero the leftmost bit set to one, with prefix 
notation. The associativity for this operator is from right to left. 

2. /: This operator is the classical not operator for binary words with prefix 
notation. The associativity for this operator is from right to left. 

3. ==: This operator is the equality operator for binary words, with infix 
notation. The associativity for this operator is from left to right. 

4. &: This operator is the classical and operator for binary words with infix 
notation. The associativity for this operator is from left to right. 

5. This operator is the classical xor operator for binary words with infix 
notation. The associativity for this operator is from left to right. 

6. =: This operator is the classical assignation operator for binary words with 
infix notation. The associativity for this operator is from right to left. 

Let us explain the use of each of these operators. 

Let 77 be a transition P system, let n be the degree of 77, let /i be the 
membrane structure of 77 and let L = {0, 1, ..., n} be the set of labels numbering 
the membranes of n, including the label 0 for the environment of the P system, 
with the restriction mentioned above for numbering the membranes of /j. The 
set of & connectivity arrays for the membrane processors system is 

{<7(z)|z€7-{0}}. (2) 

One more argument has to be incorporated to if) in order to describe the 
processing time: ip(i)(t) is the connectivity array of processor i at the evolution 
step t. 

Now let us consider what would happen during a computation in transition 
P systems. 

As it has been said before, a membrane processor changes its <7 connec- 
tivity array when it changes its operating mode from active to passive. This 
change produces changes in other membrane processors -its father and its chil- 
dren. Therefore, changes in >7’s connectivity arrays are leaded by the membrane 
processors that change their operating mode. 

Consequently, if during the transition of 77 from configuration Ct to Ct + i the 
membrane i is dissolved, then its associated membrane processor has to change 
its operating mode from active to passive, and then the following changes have 
to be produced: 

1. Membrane processor i changes its <7 (7) connectivity according to the follow- 
ing equation: 

+ 1) = ip(i)(t) A Id-Proc(i)). (3) 

2. Its father membrane processor j changes its + 1) connectivity array 

according to the following equation: 

1) = A «< ip{i){t). (4) 

Moreover, if the membrane processor j is also operating in the passive mode, 
then it has to transmit to its father membrane processor l the <7(7) (t) con- 
nectivity array in order to rearrange >7(7) (7 + 1) according to equation (4). 
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This process has to be done in a recursive manner while the father processor 
that changes its !^(Z)(£ + 1) is operating in passive mode. 

3. Its children membrane processors change their <P(k)(t + 1) according to the 
following equation: 

Vfc e L | <£■(*) (fc)(0) = 1, 

>f r (fc)(£ + 1) = <<< <f r (fc)(£) A ('?'(£)(£)& Mask-Proc(i)). (5) 

Moreover, if some of the children membrane processor fc’s are operating in the 
passive mode, they have to transmit to their children membrane processors 
Z’s the ask-Proc(i) connectivity array in order to rearrange their 

own + 1) according to equation (5). This process has to be done in a 

recursive manner while the son processors that change their \P(l)(t+ 1) are 
operating in the passive mode. 

The rest of membrane processors in the system do not change their connec- 
tivity arrays, hence keep the same \P( i ) for the next transition, that is: 

ip{i){t)=ip{i){t+l). (6) 

For all membrane processor i, the meaning of at any time is the 

following: 

— If i is in the active operating mode, then indicates the set of mem- 

brane processor labels to which i can send multisets of objects. 

— If i is in the passive operating mode, then ! V{i)(t) indicates the set of mem- 
brane processor labels to which i serves as a bridge for multisets of objects 
communication (see Fig. 1). 

Now, we shall look at another closely related issue, namely the communi- 
cation process performed by membrane processors when they have finished the 
computation of the complete multiset of evolution rules. This communication 
process also involves the connectivity array. 

Communication among membrane processors have two different ways: ascen- 
dent or bottom up, and descendent or top down. The ascendent communication 
raises no problem because there is only one membrane processor to which we 
have to communicate: the father membrane processor. However, the descendent 
communication is more complicated to implement because it could have more 
than one membrane processor to which we have to send multisets of objects. In 
order to solve this problem the cr,j binary word is defined. 

The words a,j are associated to the membrane processor i and they are 
calculated for all evolution steps for all membrane processors in the system 
associated with non elementary membranes as follows: 

V i G L - {0} 

V j G {p G L | j > i A if>(i)(p)(0) = 1}, 

c Tij = \Mask-Proc(i) & ( Id-Proc(i ) A l P(i)(t)) & <f r ( j)(t ). 



( 7 ) 
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The meaning of <jj j(k) = 1 is that the multiset of objects going from mem- 
brane processor i to membrane processor k has to be sent through the membrane 
processor j. 

We illustrate these operations with an example. Let us consider the transition 
P system 77 of degree 6 shown in Fig. 1. This figure also shows the architecture 
of membrane processors for implementing the transition P system 77 and the 
operating mode for each of the two evolution steps for the membrane processors 
associated to each membrane of 77. 






• Membrane processor in active mode 
o Membrane processor in passive mode 



Fig. 1. Membrane structure for 77 and membrane processors evolution 



The binary words that identify every processor and the initial connectivity 
array for membrane processors are the following: 

Id.Proc(l) = (0,1, 0,0, 0,0,0) 

MaskJProc{ 1) = (1, 0, 0, 0, 0, 0, 0) ip( 1)(0) = (1, 1, 1, 0, 0, 0, 0) 
Id-Proc(2) = (0,0, 1,0, 0,0,0) 

Mask-Proc{ 2) = (1, 1, 0, 0, 0, 0, 0) V’(2)(0) = (0, 1, 1, 1, 1, 0, 0) 

IdJProc{3) = (0,0, 0,1, 0,0,0) 

MaskJProc{ 3) = (1, 1, 1, 0, 0, 0, 0) V(3)(0) = (0, 0, 1, 1, 0, 1, 0) 

Id-Proc{&) = (0,0, 0,0, 1,0,0) 
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Mask-Proc{ 4) = (1, 1, 1, 1, 0, 0, 0) ^(4)( 0) = (0, 0, 1, 0, 1, 0, 1) 

Id-Proc{ 5) = (0,0, 0,0, 0,1,0) 

MaskJProc{ 5) = (1, 1, 1, 1, 1,0,0) ^(5)(0) = (0, 0, 0, 1, 0, 1, 0) 

IdJ>roc( 6) = (0,0, 0,0, 0,0,1) 

Mask-Proc{ 6) = (1, 1, 1, 1, 1, 1, 0) ^(6)(0) = (0, 0, 0, 0, 1, 0, 1) 

The transition of 77 from Cq to C\ makes disappear membrane 2 from n 
(the membrane structure of 77). Therefore the membrane processor 2 changes 
its operating mode from active to passive and its connectivity array <7(2) (1) has 
to be changed according to equation (3). Moreover, its father membrane proces- 
sor 1 has to calculate its i7(l)(l) according to equation (4). Finally, its children 
membrane processors 3 and 4 have to calculate their !7(3)(1) and *7(4) (1), re- 
spectively, according to equation (5). The rest of membrane processors do not 
change their S' connectivity arrays. For the next transition of 77, the information 
about the system membrane structure contained in membrane processors is the 
following: 

Id_Proc(l) = (0,1, 0,0, 0,0,0) ^(1)(0) = (1, 1, 1, 0, 0, 0, 0) i/>(l)(l) = (1, 1, 0, 1, 1, 0, 0) 

Id-Proc(2) = (0, 0, 1, 0, 0, 0, 0) ip( 2)(0) = (0, 1, 1, 1, 1, 0, 0) ^»(2)(1) = (0, 1, 0, 1, 1, 0, 0) 

Id-Proc(3) = (0, 0, 0, 1, 0, 0, 0) ip( 3)(0) = (0, 0, 1, 1, 0, 1, 0) = (0, 1, 0, 1, 0, 1, 0) 

Id_Proc{P) = (0,0, 0,0, 1,0,0) V*(4)(0) =(0, 0, 1, 0, 1, 0, 1) ^(4)(1) = (0, 1,0,0, 1, 0, 1) 
Id_Proc(5) = (0, 0, 0, 0, 0, 1, 0) V>(5)(0) = (0, 0, 0, 1, 0, 1, 0) ^(5)(1) = (0, 0, 0, 1, 0, 1, 0) 

Id-Proc(6) = (0,0, 0,0, 0,0,1) ^(6)(0) = (0,0, 0,0, 1,0,1) i/>(6)(l) = (0, 0, 0, 0, 1, 0, 1) 

As one can see, the connectivity array for membrane processor 1 (t/>(l)(l)) 
has been changed in order to rearrange the connections from it to membrane pro- 
cessors 3 and 4 and removing the connection with membrane processor 2. While 
the connectivity array of membrane processor 2 (^>(2)(1)) has been changed in 
order to indicate that this processor is in the passive operating mode, in fact, 
the information contained in *7(2) (1) is now related to the membrane processors 
to which processor 2 can send information through itself. 

All membrane processors in the system have to calculate their respective cr’s 
before starting the next evolution step, and this is necessary in order to update 
the system membrane structure information. In particular, for our example, at 
the end of each evolution step, the following computations have to be done for 
each of the membrane processors in order to obtain the membrane processor 
through which we send objects to their current children processors: 

— Membrane processor 1: a 12 

— Membrane processor 2: 023 and cr 2 4 

— Membrane processor 3: cr 35 

— Membrane processor 4: 

The computation of cr’s is performed according to equation (7). 

The next table shows the resulting cr’s for our transition P system 77 at the 
beginning of the transition from Cq to Cp 
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cr’s for the P system transition from Cq to C\,t = 0 


From 

(0 


Through 

(j) 


Gij =\Mask-Proc(i)&z(Id-Proc(i) A 'P(i)(t))!k'P(j)(t) 


1 


2 


<7 12 = (0, 1 , 1 , 1 , 1 , 1 , 1)&(1, 0, 1, 0, 0, 0, 0)&(0, 1, 1, 1, 1, 0, 0) = 
(0,0, 1,0, 0,0,0) 


2 


3 

4 


(723 = (0, 0, 1, 1, 1, 1, 1)&(0, 1, 0, 1, 1, 0, 0)&(0, 0, 1, 1, 0, 1, 0) = 
(0,0, 0,1, 0,0,0) 

<724 = (0, 0 , 1, 1, 1, 1, 1)&(0, 1, 0, 1, 1, 0, 0)&(0, 0, 1,0, 1, 0, 1) = 
(0,0, 0,0, 1,0,0) 


3 


5 


(735 = (0, 0 , 0, 1, 1, 1, 1)&(0, 0, 1, 0, 0, 1, 0)&(0, 0, 0, 1, 0, 1, 0) = 
(0,0, 0,0, 0,1,0) 


4 


6 


<7 46 = (0, 0, 0, 0, 1, 1, 1)&(0, 0, 1, 0, 0, 0, 1)&(0, 0, 0, 0, 1, 0, 1) = 
(0,0, 0,0, 0,0,1) 



In this case, the initial system configuration is kept by the cr’s. Objects 
communication is accomplished through the initially established communica- 
tion channel in the system configuration. Once the transition is finished, the 
membrane configuration is changed. Membrane 2 is dissolved, and then mem- 
brane processor 2 starts the rearranging process of 'P's connectivity arrays in 
order to establish the appropriate communication channels for evolution step 1. 
This evolution step starts with the C\ configuration and ends with configuration 
C 2 . The computation of !F(i)(l) have been done above, therefore, only cr’s are 
represented in the following table. 



cr’s for the P system transition from C\ to C 2 ,t = 1 


From 

(0 


Through 

(j) 


Gij =\Mask-Proc(i)&z(Id-Proc(i) A'P(i)(t))lk'P(j)(t) 


1 


2 


<7 12 = (0, 1, 1, 1, 1, 1, 1)&(1, 0, 0, 1, 1, 0, 0)&(0, 1, 0, 1, 1, 0, 0) = 
(0,0, 0,1, 1,0,0) 


2 


3 

4 


(723 = (0, 0 , 1, 1, 1, 1, 1)&(0, 1, 1, 1, 1, 0, 0)&(0, 1, 0, 1, 0, 1, 0) = 
(0,0, 0,1, 0,0,0) 

(724 = (0, 0, 1, 1, 1, 1, 1)&(0, 1, 1, 1, 1,0, 0)&(0, 1,0, 0, 1, 0, 1) = 
(0,0, 0,0, 1,0,0) 


3 


5 


(735 = (0, 0 , 0, 1, 1, 1, 1)&(0, 1, 0, 0, 0, 1, 0)&(0, 0, 0, 1, 0, 1, 0) = 
(0,0, 0,0, 0,1,0) 


4 


6 


cr 46 = (0, 0, 0, 0, 1, 1, 1)&(0, 1, 0, 0, 0, 0, 1)&(0, 0, 0, 0, 1, 0, 1) = 
(0,0, 0,0, 0,0,1) 



As one can see in the table, membrane processor 1 can send objects to mem- 
brane processor 3 and 4 through membrane processor 2 (<Ti2(3) = cr 12 (4) = 1). 
Moreover, when membrane processor 2 receives objects coming from membrane 
processor 1, it has to decide through which of its children membrane processor 
has to send such objects. In the case when the objects go to membrane proces- 
sor 3 it will send the objects through itself, because it is in the active operating 
mode (< 723 ( 3 ) = 1). In the case when the objects go to membrane processor 4 
it will also send the objects through itself, because it is in the active operating 
mode ((724(4) = 1). 
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Now, one more transition of 77 is accomplished in our example, and the P sys- 
tem goes from C\ to Ci dissolving membranes 3 and 4. Once again, the *7(3) (2) 
and !7(4)(2) connectivity arrays of membrane processors have to be changed us- 
ing equation (3), while membrane processors 2 and 1 change their ip( 2) (2) and 
^(l) (2) connectivity arrays using equation (4) and membrane processors 5 and 
6 do it using equation (5). The resulting system I 1 connectivity arrays for this 
evolution step ( t = 2) are the following: 

Id_Proc( 1) = (0, 1, 0, 0, 0, 0, 0), (1)(0) = (1,1, 1, 0, 0, 0, 0), V>( 1)(2) = (1, 1, 0, 0, 0, 1, 1), 

Id_Proc{ 2) = (0, 0, 1, 0, 0, 0, 0), ^(2)(0) = (0, 1,1,1, 1, 0, 0), V>(2)(2) = (0, 1, 0, 0, 0, 1, 1), 

Id_Proc( 3) = (0,0, 0,1, 0,0,0), ip{ 3 ){ 0 ) = (0, 0, 1, 1, 0, 1, 0), V>(3)(2) = (0, 1, 0, 0, 0, 1, 0), 

Id_Proc( 4) = (0, 0, 0, 0, 1, 0, 0), </>(4)(0) = (0, 0, 1, 0, 1, 0, 1), V>(4)(2) = (0, 1, 0, 0, 0, 0, 1), 

Id.Proc( 5) = (0,0, 0,0, 0,1,0), ip( 5 ){ 0 ) = (0, 0, 0, 1, 0, 1, 0), ip( 5 ){ 2 ) = (0, 1, 0, 1, 0, 1, 0), 
Id_Proc( 6 ) = (0,0, 0,0, 0,0,1), ^(6)(0) = (0, 0,0,0, 1,0,1), -0(6)(2) = (0, 1,0,0, 1,0, 1). 

After this transition, the membrane system has changed its membrane struc- 
ture once again, hence the membrane processors have to recompute their 0’s 
using equation (7). The next table shows these values for membrane processors 
at t = 2. 



0’s for the P system transition from C2 to C3,t = 2 


From 

(0 


Through 

(j) 


c Tij =\Mask-Proc(i)&i(Id-Proc(i) A <7(i)(t))&!7(j)(t) 


1 


2 


0 12 = (0, 1, 1, 1, 1, 1, 1)&(1, 0, 0, 0, 0, 1, 1)&(0, 1, 0, 0, 0, 1, 1) = 
(0,0, 0,0, 0,1,1) 


2 


3 

4 


023 = (0, 0, 1, 1, 1, 1, 1)&(0, 1, 1, 0, 0, 1, 1)&(0, 1, 0, 0, 0, 1, 0) = 
(0,0, 0,0, 0,1,0) 

024 = (0, 0, 1, 1, 1, 1, 1)&(0, 1, 1, 0, 0, 1, 1)&(0, 1, 0, 0, 0, 0, 1) = 
(0,0, 0,0, 0,0,1) 


3 


5 


035 = (0, 0, 0, 1, 1, 1, 1)&(0, 1, 0, 1, 0, 1, 0)&(0, 1, 0, 0, 0, 1, 0) = 
(0,0, 0,0, 0,1,0) 


4 


6 


0 46 = (0, 0, 0, 0, 1, 1, 1)&(0, 1, 0, 0, 1, 0, 1)&(0, 1,0,0, 0, 0, 1) = 
(0,0, 0,0, 0,0,1) 



After the last transition, membrane processor 1 can send objects to mem- 
brane processors 5 and 6. Once again, this has to be done through membrane 
processor 2 (012(5) = 012(6) = 1). When membrane processor 2 receives ob- 
jects from membrane processor 1, it has to decide through which of its children 
membrane processor has to send the objects. In the case when the objects go to 
membrane processor 5 it will send through membrane processor 3 (023(5) = 1), 
otherwise it will send through membrane processor 4 (024(6) = 1). When mem- 
brane processors 3 or 4 receive objects from membrane processor 2 -they are 
in the passive operating mode- they have to decide through which of its chil- 
dren membrane processor will send such objects. In both cases, they only have 
one child membrane processor, which is in the active operating mode. There- 
fore they will send the objects to their corresponding child membrane processor 
(0-35(5) = 1 and 046(6) = 1). 
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The membrane system is now ready to start a new transition (although the 
example does not discuss it). 

4 Conclusions 

This paper illustrates the membrane processors operational mode based on using 
connectivity arrays and the way to update the system membrane structure. It 
also illustrates the work of membrane processors in order to solve the commu- 
nication problems caused by membrane dissolution in transition P systems. 

It has also been shown that additional information about the membrane 
structure is needed in order to obtain the desired functionality in membrane 
processors for solving communication problems. Binary words and the defined 
operations on them constitute a well defined data structure that provides the 
desired functionality to the membrane processor systems. 
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Abstract. New classes of P automata are introduced corresponding to 
the basic classes of languages in the Chomsky hierarchy. Unlike the pre- 
viously defined P automata, active P automata are computing with the 
structure of the membrane systems, using operations like membrane cre- 
ation, division and dissolution. The model is applied to the parsing of 
(natural language) sentences into dependency trees. 



1 Introduction 

P automata are a special type of membrane systems that work as accepting de- 
vices. Different classes of P automata have been introduced and studied recently, 
in papers as [3,5,6,12,1,4,11]. 

In [3,4], one defines P automata using communication rules and final config- 
urations. In [1], there are used both evolution and communication rules. In [11], 
again evolution and communication rules are used, but the input symbols are 
not introduced in the outer membrane as in the first two cases but they exist on 
a tape which is transfered between membranes during the computation. 

In all these approaches, the P automaton preserves the same structure of 
membranes during the whole computation. No rules of membrane creation, divi- 
sion or dissolution are involved. With active P automata, we propose a parallel 
computation model with a flexible structure that evolves in time according to 
the set of rules. 

Informally, in active P automata the accepting computation starts with one 
membrane, which contains the string that has to be recognized and possibly some 
other information. The computation follows according to the input string and 
during the evolution membranes can be created or dissolved. The computation 
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